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METHOD AND SYSTEM FOR REPRESENTATION OF 
GEOGRAPHICAL FEATURES IN A COMPUTER-BASED SYSTEM 

TECHNICAL FIELD OF THE INVENTION 
The present invention relates generally to cartographic databases, and in particular, to the 
representation of geographic features using polynomial splines. 

BACKGROUND OF THE INVENTION 
Computer-based systems that display graphical maps of geographic areas use 
cartographic databases that contain data that represent the geographic features in the areas. The 
geographic features represented by data in cartographic databases may include roads, 
intersections, points-of-interest, lakes, railroad tracks, buildings, airports, stadiums, parks, 
mountain ranges, bridges, and so on. To permit these features to be displayed, a cartographic 
database includes data that represent the geometry (e.g., size, shape and position) and character 
of these various geographic features. 

In order to represent roads (as well as other geographic features on a computer), it is 
necessary to use some form of discretization. In some cartographic databases, each road is 
represented in terms of one or more discrete road segments. Each road segment is represented 
by a separate data entity contained in the cartographic database. Thus, each road segment data 
entity represents a portion of a roadway in a particular geographic region. 

In some cartographic databases, the geometry of each represented geographic feature is 
represented using one or more discrete data points wherein each data point consists of coordinate 
pairs (longitude, latitude) or coordinate triples (longitude, latitude, altitude). A point feature, 
such as a point of interest or intersection of roads, is represented using a single data point, while 
continuous features, such as road segments, lake boundaries, and so on, are represented as a 
sequence of data points. The end points of the sequence of data points used to represent a 
continuous geographic feature are referred to as nodes, while data points between the end points 
are referred to as shape points. Using nodes and shape points, the true shape of a continuous 
geographic feature can be approximated. The shape points and nodes can be stored in a shape 
point database, which may be part of a larger cartographic database. 



1 



To display the geographic features on a computer display device, piecewise linear 
interpolation is used between the shape points to draw straight lines connecting the points. This 
is illustrated in FIG. 1, which depicts a geographic feature (solid line) approximated by linearly 
5 interpolated lines (dashed lines) connecting shape points SP(0) - SP(5) and the nodes. With this 
representation scheme, the spacing between shape points is selected such that the error between 
the data point representation and the actual geographic feature lies within some acceptable 
bound. 

While using data points and linear interpolation to represent the shape of geographic 
10 features has the advantage that it is conceptually relatively simple to work with in database 
development, compilation, access and display software, there is opportunity for improvement. 
For example, with linear interpolation, to achieve a high level of accuracy and detail, it is 
sometimes necessary to store a relatively large number of data points for a given feature. This 
increases the storage requirements of the cartographic database, as well as the database read 
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Accordingly, there is a need for a method and system of representing and storing 
geographic features, such as road geometry, that reduces the storage requirements while 
providing a high level of accuracy in representing the features. Further, it is desired that the 
increases in resolution and accuracy in the cartographic database do not negatively impact the 
database read time. 

H Another way to represent the geometry of continuous geographic features in a 

cartographic database is disclosed in U.S. Patent No. 6,029,173. The '173 patent describes a 
navigation system having a cartographic database that stores data points and Bezier control 
points that represent road geometry. Although the approach disclosed in the '173 patent provides 
25 an advantage over prior methods for representing geographic features, there continues to be room 
for improvement. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 depicts an exemplary geographic feature represented by shape points and nodes 
30 with piecewise linear interpolation. 

FIG. 2 depicts the geographic feature of FIG. 1 represented by a polynomial spline. 
FIG. 3 depicts a parametric polynomial spline. 
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FIG. 4 is a block diagram illustrating a system for converting a shape point database into 
a spline database according to an embodiment of the present invention. 

FIG. 5 is a flow chart illustrating a method of operating the processor shown in FIG. 4. 

FIG. 6 depicts an exemplary geographic feature having a straight section represented by 
shape points and curved sections represented by spline control points. 

FIG. 7 shows a geographic feature data entity for storing shape points and control points 
representing a geographic feature having both curved and straight sections. 

FIG. 8 is a block diagram illustrating a computer-based system capable of displaying 
geographic data, such as roadways represented by polynomial splines in accordance with a 
further embodiment of the present invention. 

FIG. 9 is a flow chart illustrating the operation of the computer-based system shown in 

FIG. 8. 

DETAILED DESCRIPTION OF THE PRESENTLY 
PREFERRED EMBODIMENT(S) 

The present invention provides a polynomial sphne based method for representing 
geographic features in cartographic databases. The use of polynomial splines reduces the amount 
of data needed to accurately represent geographic features in cartographic databases. This results 
in smaller database sizes and shorter database read times, while maintaining or enhancing the 
accuracy with which geographic features are depicted in the database. 

Turning now to the drawings, and in particular to FIG. 2, there is illustrated a polynomial 
spline representation of the geographic feature shown in FIG. 1 . In FIG. 2, the spline is shown 
as a dashed line, while the actual geometry of the geographic feature is shown as a solid line. 

A polynomial spline consists of piecewise polynomial functions which are joined end to 
end with some continuity condition imposed at the join points. A polynomial spline is 
determined completely by a set of control points CP(0) - CP(N). and in some cases, additional 
parameters (e.g., knots and/or control point weights). 

Splines can be broadly classified as interpolating and noninterpolating. For both types of 
splines, the control points can be computed so that a spline accurately represents a corresponding 
geographic feature. 



A polynomial spline refers to any composite curve formed with piecewise polynomial 
functions representing curve sections and satisfying specified conditions at the boundaries of the 
sections. There are many types of spline curves. 

FIG. 3 depicts a planar parametric spline curve. The intermediate points k(0) - k(n) 
w^here the spline sections meet are known as knots. In general, to ensure that the mth derivatives 
of the spline are continuous at the knots, a spline of order m+1 or greater is used. When the 
knot spacing is uniform, the spline is said to be uniform, and conversely, when the knot spacing 
is non-uniform, the spline is termed as a non-uniform spline. 

I. Database Conversion 

FIG. 4 is a block diagram of a system 10 for converting a database that represents feature 
geometry using shape points and nodes into a database that represents feature geometry using 
spline control points. The system includes a processor 12 in communication with a shape point 
database (DB)14 and a spline database 16. The shape point database 14 contains shape points 
and nodes. A user interface 1 1 permits users to configure conversion parameters used by the 
processor 12. The shape point and spline databases 14,16 can be stored on any suitable 
computer-usable storage medium, such as one or more CD ROMs, hard drives, DVDs, diskettes, 
solid-state memories such as RAMs, ROMs, or EEPROMs, or the like. As will be discussed in 
greater detail below, the processor 12 processes data stored in the shape point database 14 to 
convert it into a plurality of spline control points that are stored along with other defining 
parameters in the spline database 16. The spline control points can provide a more efficient 
representation of map data in terms of storage requirements. 

The shape point database 14 stores a plurality of shape points and nodes representing one 
or more geographic features that can be displayed by a computer-based system, such as a 
navigation system, as part of a digital map. Shape points and nodes represent geographic 
locations within a map. In the embodiments disclosed herein, the shape points and nodes can 
consist of coordinate triples (longitude, latitude, altitude) or coordinate pairs (longitude, latitude), 
which represents specific locations along geographic features. Nodes are special cases of shape 
points, representing, for example, intersections between road segments and end points of 
geographic features. The shape points can be part of a cartographic database used for map 
display. The cartographic database itself may be part of a larger geographic database, which 



may include other kinds of databases used for functions such as route calculation, maneuver 
generation, and other navigational functions. 

Data stored in the shape point database 14 can be organized on a per geographic feature 
5 basis, i.e., each geographic feature is represented as sequence of nodes and shape points. 
Although the database 14 can store a number of geographic feature entities 17, three entities 
SEGo, SEGi, SEG2 are shown for illustration purposes. A pair of nodes 23 and zero or more 
shape point entries 25 can be associated with each entity. For entities representing straight or 
substantially straight geographic features, the nodes 23 are sufficient to accurately represent the 
10 feature; and therefore, for straight geographic features, the nodes alone can be stored without 
shape points. 

In the example shown, the entities 17 correspond to geographic features (e.g., roadways) 
displayable on a map 15. A first geographic feature SEGo is represented by shape points SPo(O) 
^5 - SPo(i), where i is an integer, and nodes ENo(O) - ENo(l); a second feature SEGi is represented 

15m by shape points SPi(O) - SPiG), where j is an integer, and nodes ENi(O) - ENi(l); and a third 

ill 

Iq feature SEG2 is represented by shape points SP2(1) - SP2(k), where k is an integer, and nodes 

fl EN2(0)-EN2(1). 

ifl The processor 12 can be any computing device for performing the steps shown in the 

L, flow chart 20 of FIG. 5, including a conventional personal computer (PC) executing a software 

2!^ program. In such an embodiment, the software program can be stored on a computer-usable 

y 

fjl media for configuring the processor 12 to perform the method of FIG. 5 as disclosed herein. The 
■ij computer-usable media can be any article of manufacture for storing computer programs, such as 
a CD ROM, hard drive, DVD, diskette, solid-state memory such as a RAM, ROM, or EEPROM, 
or the like. Alternatively, the processor 12 can be a custom hardware device, such as one or 
25 more application specific integrated circuits (ASIC) designed to perform according to the method 
shown in FIG. 5. 

Generally, for each geographic feature entity in the shape point database 14, the 
processor 12 generates a corresponding set of control points and other parameters 19 that define 
the spline; these parameters are stored in the spline database 16. These control points and other 
30 parameters can be used to compute the spline function, which can represent the geographic 
feature in a displayable map 21. 
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In the example shown, the first geographic feature SEGo is represented by control points 
CPo(O) - CPo(l), where 1 is an integer; the second geographic feature SEGi is represented by 
control points CPi(O) - CPi(m), where m is an integer; and the third geographic feature SEG2 is 
5 represented by shape points CP2(0) - CP2(n), where n is an integer. 

FIG. 5 shows a flow chart 20 of a method of operating the processor 12 in accordance 
with an embodiment of the present invention. In step 22, the processor 12 receives shape points 
and nodes representing one or more geographic features from the shape point database 14. The 
processor can include a memory (not shown) for buffering the shape points. 
10 Next, in step 24, the number of spline control points to be generated for each segment is 

determined. Generally, the number of control points per segment is set to be less than the 
number of shape points per segment. This results in smaller database storage requirements. The 
number of control points per segment can be set to a predetermined value or fraction of the 
number of shape points, or alternatively, the number can be selected by a user or be dynamically 

154 computed on a per segment basis depending on the desired level of representational accuracy. 

fti 

.g The user interface 1 1 can be used to select the number of control points either globally, for all 

'^^ geographic features being processed, or individually for each of the geographic features. Other 

W 

\n parameters which are dependent on the type of spline being computed are also determined, either 
globally or on a per geographic feature basis. 

2W In step 26, a least squares approximating polynomial spline is fitted to each of the 

0 

in geographic features. The polynomial spline can be any type of spline, for example, a uniform or 
nonuniform nonrational B-spline, NURBS, or a uniform or nonuniform CatmuU-Rom spHne. 
The application of each of these types of splines is described in further detail below. 

The spline is fit to the shape points and nodes representing each geographic feature by 
25 applying a least squares approximation to the shape points and other parameters that define the 
geographic feature to generate the selected number of spline control points for each segment. 

As will be discussed in further detail below, regularization techniques can be employed to 
make the least squares approximation better conditioned. In addition, "bearing values" and 
selective weighting can be included in the least squares approximation to improve the continuity 
30 of spline curves between geographic features that meet at a point. 

In step 28, the spline control points for the geographic features are then stored in the 
spline database 16 (step 28). 



6 



II. Generalized Approach to Computing Least Squares Approximating Polynomial Splines 
Although the present invention is not limited to any particular method of performing a 

least squares approximation, this section describes a least squares methodology in accordance 
5 with one embodiment of the invention. The least squares methodology described in this section 

can be used with any type of spline function, including uniform B-splines, non-uniform B- 

splines, or CatmuU-Rom splines. 

To fully appreciate the least squares method described herein, it is appropriate to first 

provide a general mathematical description of a polynomial spline and its representation in 
10 parametric form. 

A polynomial spline space curve of any type can be represented in parametric form as a 

set of polynomial functions of an independent parameter / as follows: 

O 

I H()=M) Pyit) pM (1) 

ri 

1^ where 

*P 1=0 
M A/-1 

III 

M PAo=l,p;N,ity, (2) 

and X , y and z represent longitude, latitude and altitude, respectively. The set of functions 
A^-(0 are called spline basis functions, and the weights P/" , P/" and P'on these basis functions 
are the control points of the spline. The basis functions vary for different types of splines. In the 

25 case of uniform splines, e.g., uniform B-spline and uniform CatmuU-Rom spline, the set of basis 
functions is obtained by shifting a single basis function over the length of the geographic feature. 
In the case of nonuniform splines, e.g., nonuniform B-splines, the set of basis functions is 
controlled by the knot locations. For NURBS, the most general form of B-splines, the shape of 
the approximating functions is controlled by the knot locations, as well as by weights. The knot 

30 locations can be selected by a user for each geographic feature via the user interface 11. 
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Each geographic feature is specified by a set of shape points and nodes, where each point 
is specified as a (longitude, latitude, altitude) coordinate triple. Accordingly, for each 
geographic feature there is a list of points (xo>>^o'^o)» {^\^y\y^\)^ {^N-\^yN-\^^N-\) where 
the Xi are the longitude values, the y. are the latitude values and the z. are the altitude values. 
The first and last coordinate triple, viz., {xQ^yQ^z^) and (^Ar-i'>^yv-P^A^-i) nodes, while 

the coordinates between these values are shape points. 

In order to obtain a polynomial representation of the shape points, it is necessary to 
parameterize the given shape points in terms of a parameter / . This is accomplished by 
assigning a value t- to each triple {x.,y.,Zi), 

Parameterization can be accomplished in a number of different ways. Parameterization 
schemes can be broadly classified as uniform and nonuniform. In the first case, the assigned 
parameter values /. are uniformly spaced, generally with tQ=0 and each succeeding value 

exceeding the previous value by some constant. 

A particularly successful nonuniform parameterization method involves assigning = ^ 

to a first shape point and incrementing the parameter value by the distance between the first 
shape point and the next shape point. This is called chord length parameterization. 

For each geographic feature, parameterization results in a set of parameter values 
ro,/p...,/;y_i corresponding respectively to the latitude values x^,x^,,„,x^^^ , longitude values 

and altitude values Zo,Zi,...,z^., . The parameter values and correspond to 
the nodes and the parameter values ,...,/;v_2 correspond to the shape points. 

To reduce the error in the least squares fit polynomial spline, "bearing values" are 
determined for the nodes of each geographic feature. The bearing values are the first derivatives 
at the nodes, i.e., at / = and at / = . These values indicate the direction in which the 

geographic feature is headed at its endpoints. 

The purpose of the least square fit is to find the optimal parametric polynomial spline 
curves P^it), Py{t) and P^{t) such that the mean squared error between (/) and the given 

longitude and bearing values, Py {t) and the given latitude and bearing values and P^ (t) and the 

given altitude and bearing values at the parameter values ^0,/^,...,/^.^ is minimal. It is preferable 
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that the nodes of each geographic feature have a larger weight associated with them than the 
shape points so that the spHne is forced closer to the nodes. This is important so that continuity 
between connected geographic features can be maintained in the spline representation. 

Since the steps for determination of P^(t)y Py{t) and P^{t) are identical, the least 

squares methodology is illustrated for P(/) , where P(0 represents P^{t), Py{t) and P^(t), This 

simplifies the notation somewhat for illustrative purposes. 

Accordingly, the methodology defines a function ^^O) follows 



r{Pit)) = I [P{tj ) - V, ) + W{P{t, )-v,Y^ W{P(t,,, ) - v^., y 

+ (^'(^o)-v'(/o))'+(P(^^_,)-v'(r^.O)' 



(3) 



where Wis the weight assigned to the nodes and P(t) is understood to represent P^(t) , Py(t) 

and P^ (0 and v is understood to represent x , >' or z . 

From Equafion (3), P(t) can be expressed in terms of its control points and the basis 
function of the spline as follows: 



Substituting Equation (4) into Equation (3), yields Equation (5) below. 

N'2fM-\ y y (m-\ \ 



(4) 



7=1 V 1=0 



\i=0 



+ 



V /=0 J \ i=0 J 



(5) 



The least squares approximation is computed by minimizing r(P(0) of Equation (5). 
However, the resulting minimization is not in general "well conditioned", i.e., in solving the 
system of equations represented by Equation (5), the system matrix may be close to singular. 
There is therefore a need to use a regularization term, which will make the problem better 
conditioned. 



Although various forms of regularization can be used in accordance with the present 
invention, one embodiment of the invention employs a specialized form of Tikhonov-Miller 
regularization, which in general involves placing a constraint on some function of P{t) and its 
derivatives. In particular, the third derivative of P{t) can be used, though this is only one 
possibility. Accordingly, the regularization modifies Equation (5) as follows: 



N-2 







2 / 






+ w 


V/=o 


) 





fM-x y (m-\ V 



Vi=o 



'o 



(6) 



N-2f M-\ 

is- 

7=1 V /=0 



Vi=o 



. /=0 



^ /=0 J V /=0 J u \ i=0 



V fM-\ 
J \ i=0 



2 r 



(7) 



dt 



where /is a constant which determines the degree of regularization added. 

As stated above, the goal is to find the P(0 that minimizes r{P{t)). Since P(0 is a 
function of the control points P^,P^,...,P^f_^ , it is therefore the goal to find the set of control 
points Po, which minimizes r{P{t)). To accomplish this, r(P(0) is differentiated 
with respect to P, , and the resulting derivative is set equal to zero, as shown in Equations (8) - 
(9). 
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^2w['f^P,N,{tj)-v}N,itj) + 2w['f^P,N,{tj)-v,\^ 

V/=o J V/=o / 



+ 2 



or 



\i=0 J Vi=0 / (o V'=0 



A^, '{t)dt 



(8) 



A/-1 N-2 . , , . ^ 

r \ 



A/-1 
/=0 



^-2 



r \N;{t)N;{t)dt 



= S^^/(0)-^^oA^,(^o)-^v^-.A^,0;v-,)-v"(^oW/'(^)-v'(^;v-.W'(^;) 



(9) 

Equation (9) represents one row of a system of equations which is solved for the unknown 
coefficients P^,P,,,..,P^_,. Differentiating r(P(0) with respect to P,,/ = 0,...,M-1, and setting 
each result equal to zero, the full system matrix is obtained. Since the goal is to determine (t) , 
P (0 and P^ (t) 5 the following systems of equations are solved of the form of Equation (6): 



A/-1 f l^-l , . . . ^ 



M-\ 



= Z^j^iih^- Wx,N,(t,)- Wx^_,N,{t,_,) - x\t,)N;itj)-x\t^.,)N;itj) 



(10) 
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/=0 



N-2 



V ^0 



= T,yjN,itj)-Wy,N^{t,)-Wy^_,N^it,,,)-y\i,)N;it^ 



and 



(11) 



M-l f N-1 \ 

Y.N,{tj)N^{tj)^WN,{t,)N,{t,) + WN,{t,_,)N,{t,.,)^N 

i=0 V J 



M-l 



= f^ZjN,{_tj)-Wz,N,{t,)-Wz^.,N,{t,.,)-z{t,)N;{tj)-z\t,.,)N;{tj) 



(12) 



Accordingly, the polynomial spline coefficients are determined by solving Equations 
(10H12). 

As described earlier, uniform or nonuniform nonrational B-splines, uniform or 
nonuniform Catmull-Rom splines, NURBS, and the like can be used to represent geographic 
features. The details specific to some of these spline types are given below. 



A. Uniform Nonrational B-Spline Representation 

Equations (10) - (12) are solved for the spline coefficients ,P^^ ,,..,P^_^ , 

and Po',f;',...,P^_i. To solve these equations, the / th 4^%rder uniform B-spline 

basis function A^.(/) is substituted into Equations (10) - (12). The uniform nonrational B-spline 

basis function N^{f) is defined as follows: 



12 



«,(')- 



/ - 3 < / < / - 2 



- 1/^ + i/' + 1/ + 1, / - 2 < / < / - 1 



2 6 



1 3 2 2 

— - r + — , 

2 3 



/-!</</ 



6 2 2 6 



0 otherwise 



(13) 



Therefore, 

W 

pi 

3 . 

m 

In 



1, i-3<t<i-2 
i-2<t<i-\ 



-3, 



3, 



1, 



0 otherwise 



(14) 



Now consider the regularization term, marked / in Equation (9) - (12). Because the 
10 third derivative of of a 4^^ order polynomial is a constant, the integral / can be to computed: 



I = |Ar;" (OTv/" {t)dt = Ja^;' (OiV/"' (0^^ + J (0^1^ (Odt + - • + j^, (o a^/ (o^^ 



(15) 
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But A^/"(0 has support only from i - 3 < / < / + 1 . Therefore, 

/= \N:"it)N,''{t)dt=o+-+ lN.''it)N;\t)dt+ iN:\t)N;\t)dt+ |<"(OA^/'"(Orf^+ jN'(t)N;"it)dt 

tQ /-3 i-2 i-\ 

i-2 i-\ i (+1 

= \N"{t)dt-3\N^\t)dt + 3\Ni\t)dt- \Ni'{t)dt 

i-3 i-2 1-1 I 

(16) 

where some of the terms may be zero, depending on the value of / . Equation (16) 
represents the regulation term used with the uniform B-spline. 

B. Uniform Catmull-Rom Spline Representation 

Like the uniform B-spline, Equations (10) and (12) are solved for the spline coefficients 
coefficients , P,\P,\,.,,P^_, and . However, the /th4^%rder 

uniform Catmull-Rom spline basis function N-{t) is used instead of the uniform B-spline basis 
function. The Catmull-Rom basis function is defined as follows: 



N,{t) = 



1 3 1 2 

2 2 



i-Z<t<i-2 



--t'+2t^+-t, i-2<t<i-\ 



2 2 

i-\<t<i 



1 , ,1 

— f +r — t, i<t<i + \ 

2 2 

0 otherwise 



(17) 



Therefore, 
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si 



u 



3, i-3<t<i-2 
-9, i-l<t<i-\ 

-3, i<t<i + \ 
0 otherwise 



The regularization term / is within a constant of the corresponding regularization term 
for the uniform B-spline described above. However, since the basis functions for this spUne are 
5 different, the system of equations to be solved will have a different system matrix and known 
vector. 

C. Nonuniform Nonrational B-Spline Representation 

The basic set of systems to be solved is Equations (10) and (12). However, the 
0 regularization term / is not necessarily required in this case. Further, the knot values can be 
1 selected by the user. 

iu 

^0 D. NURBS 



Again, the basic set of systems to be solved is Equations (10) and (12). As for 
nonuniform nonrational B-splines, the regularization term/ may not be required in this case. 



1^ Further, the knot values as well as control point weights can be selected by the user, thus 



n affording the user greater control over the shape of the basis function and the spline. 



III. Polynomial Spline with Alternative Representation for Straight Road Sections 

FIG. 6 depicts an exemplary geographic feature having a straight section represented by 

20 shape points and curved sections represented by a polynomial spline. For geographic features 
having long straight sections, the straight sections can be represented using the conventional 
linear interpolation method. The shape point coordinates of the endpoints of the straight section, 
the control points corresponding to the curved sections and other defining parameters can be 
stored in a single geographic feature entity in the spline database 16. 

25 The curved sections can be represented using a best fit polynomial spline. Continuity 

is ensured at the point where the straight section and curved sections of the segment meet by 
incorporating this information into the least squares approximation of the curved sections of the 
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geographic feature. Further, the degree to which this continuity is enforced can be controlled in 
the least squares approximation by the configurable weight factor. 

The methodology to accomplish this is similar to that described above for polynomial 
splines, except that each geographic feature is split into sections. At points where a spline 
section meets a straight line section, the bearing values, i.e., first derivatives, of the two sections 
are incorporated into the least squares approximation. In other words, one of the bearing values 
in Equations (10) and (12) is replaced by the bearing value of the straight line segment. In this 
manner, straight sections of a geographic feature can be represented using conventional shape 
points and linear interpolation, while curved sections of the same feature can be represented 
using polynomial spline control points. 

FIG. 7 shows an exemplary geographic feature data entity 40 for storing shape points 44, 
control points 43, and other defining parameters representing a geographic feature having both 
curved and straight sections. The data entity 40 can be stored in the spline database 16. 

Included in the entity 40 is an alternative representation (AR) flag 42. The AR flag 42 
can indicate to a display system reading the entity 40 that the entity 40 is represented using both 
control points and shape points. Accordingly, display systems can appropriately perform either 
linear interpolation or spline calculation using the correct entries in the entity 40. 

IV. Display System for Polynomial Splines and ADASs 

FIG. 8 is a block diagram illustrating a computer-based system 50 capable of displaying 
geographic data, such as roadways represented by polynomial splines. The system 50 includes a 
computer 52, a spline database (DB) 54, and an output display device 56. 

The computer 52 includes a processor 58 and a memory 60. The processor 58 can 
execute a computer software program stored in the memory 60. The software program can direct 
the processor 58 to perform the steps shown in flow chart 80 of FIG. 9. 

The flow chart 80 illustrates the operation of the computer-based system 50 shown in 
FIG. 8. In step 82, the processor 58 retrieves spline control points and other defining parameters 
from the spline database 54. Next, in step 84, the processor 58 calculates the polynomial spline 
representing a geographic feature. Through a conventional display interface card, the processor 
58 then passes the spline information to the output display device 56, which then displays the 
spline representing the geographic feature to the user (step 86). 
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The spline database 54 can include the same data included in the spline database 16 
shown in FIG. 4. 

The polynomial spline calculated by the processor 58 can be any type of spline, including 
a uniform or nonuniform nonrational B-spline, NURBS, and a Catmull-Rom spline. In addition, 
the computer 52 can be configured to display geographic feature using the alternate 
representation technique described above for geographic feature having both straight sections 
and curved sections. 

The computer-based system 50 of FIG. 8 can be used in navigation systems included with 
on-road vehicles, such as automobiles, trucks, buses, and so on, as well as advanced driver 
assistance systems (ADASs) provided in on-road vehicles. 

Advanced driver assistance systems are capable of improving the safety, comfort, 
efficiency, and overall satisfaction of driving. Examples of advanced driver assistance systems 
include adaptive headlight aiming, adaptive cruise control, and adaptive shift control. Adaptive 
headlight aiming adjusts the vehicle's headlights, i.e., width, rotational angle, elevation angle, 
and brightness, based on the curvature of the road ahead of the vehicle, tilt, elevation change, and 
other factors. Adaptive cruise control allows the vehicle to maintain and/or resume a set speed 
or safe following distance from other vehicles at slower than the set speed based on data about 
vehicle speed, nearby vehicles and other obstructions, type of road traveled (motorway vs. local 
road), road curvature, tilt, elevation, and other factors. Adaptive shift control adjusts the gearing 
and shifting of automatic transmissions based on sensor data about vehicle speed, engine speed, 
road curvature, tilt, elevation, and other factors. There are other advanced driver assistance 
systems in addition to these. 

These advanced driver assistance systems use a variety of sensor mechanisms in the 
vehicle to determine the current state of the vehicle and the current state of the roadway in front 
of the vehicle. These sensor mechanisms may include radar and vision-oriented sensors, such as 
cameras. Although radar and vision-oriented sensors are important components of ADASs, these 
components do not always provide a completely accurate representation of driving conditions. 
For example, the range and/or accuracy of radar or vision-oriented sensors can be affected by 
certain environmental conditions, such as fog, heavy rain or snow, or snow-covered roads. 
Moreover, radar and vision-oriented systems do not reliably detect certain usefiil road attributes. 
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such as speed limits, traffic signs, bridge crossings, etc. Further, radar and vision-oriented 
sensors cannot "see" around comers or other obstructions and therefore may be hmited under 
such circumstances. 

One way to improve radar and vision-oriented systems is to use digital map data, such as 
road geometry information, as an additional component in ADASs. Digital map data can be used 
in ADASs to provide information about the road ahead. Digital map data are not affected by 
environmental conditions, such as fog, rain or snow. In addition, digital map data can provide 
useful information that cannot reliably be provided by vision-oriented systems, such as speed 
limits, traffic and lane restrictions, or the like. 

Further, road geometry information included in the digital map data can be used to 
determine the road ahead of the vehicle even around comers or beyond obstmctions. 
Accordingly, digital map data, and in particular, road geometry information, can be a useful 
addition in ADASs. 

Representing road geometry using polynomial splines, as described herein, provides an 
accurate and compact way of representing such information in databases suitable for use in 
ADASs. These databases can improve the overall performance of ADASs by increasing the 
accuracy of the digitized road geometry used by the ADASs, while reducing the storage space 
required for the digital map information. 

While specific embodiments of the present invention have been shown and described, it 
will be apparent to those skilled in the art that the disclosed invention may be modified in 
numerous ways and may assume many embodiments other than those specifically set out and 
described above. Accordingly, the scope of the invention is indicated in the appended claims, 
and all changes that come within the meaning and range of equivalents are intended to be 
embraced therein. 
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